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ABSTRACT 

Clusters of galaxies have been used extensively to determine cosmological parameters. A major difficulty in 
making best use of Sunyaev-Zel'dovich (SZ) and X-ray observations of clusters for cosmology is that using 
X-ray observations it is difficult to measure the temperature distribution and therefore determine the density 
distribution in individual clusters of galaxies out to the virial radius. Observations with the new generation of 
SZ instruments are a promising alternative approach. We use clusters of galaxies drawn from high-resolution 
adaptive mesh refinement (AMR) cosmological simulations to study how well we should be able to constrain 
the large-scale distribution of the intra-cluster gas (ICG) in individual massive relaxed clusters using AMiBA 
in its configuration with 13 1.2-m diameter dishes (AMiBA 13) along with X-ray observations. We show 
that non-isothermal j3 models provide a good description of the ICG in our simulated relaxed clusters. We 
use simulated X-ray observations to estimate the quality of constraints on the distribution of gas density, and 
simulated SZ visibilities (AMiBA 13 observations) for constraints on the large-scale temperature distribution 
of the ICG. We find that AMiBA13 visibilities should constrain the scale radius of the temperature distribution 
to about 50% accuracy. We conclude that the upgraded AMiBA, AMiBA 13, should be a powerful instrument 
to constrain the large-scale distribution of the ICG. 

Subject headings: galaxies: clusters: general 



1. INTRODUCTION 

According to our standard structure formation scenarios 
based on the dark matter (DM) models, clusters of galax- 
ies, the most massive virialized objects in the Universe, form 
from the largest positive density fluctuations, thus their for- 
mation and evolution are sensitive to the underlying cosmo- 
logical model. Taking advantage of this feature, clusters have 
been used extensively to determine cosmological parameters 
(e.g., Henry 2000; Allen et al. 2004; Ettori 2004; Vikhlinin 
et al. 2008; for recent reviews see Voit 2005 and Borgani 
2006). Prospects of determining cosmological parameters us- 
ing much larger samples of clusters of galaxies from next 
generation surveys were discussed in detail by e.g., Haiman, 
Mohr & Holder (2001); Holder, Haiman & Mohr (2001); 
Molnar, Birkinshaw & Mushotzkyl f2002 ) and IMolnar et al. 
(2004jT~ 

While theory predicts the mass function of clusters of 
galaxies and the distribution of mass, gas density and temper- 
ature within individual clusters, observations directly measure 
the X-ray luminosity and intra-cluster gas temperature func- 
tions, luminosity functions based on the Sunyaev-Zel'dovich 
(SZ) effect, and the projected distribution of X-ray emissiv- 
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ity and electron pressure. To connect theory and observation 
it is crucial to understand the physics of clusters out to their 
virial radii and beyond. The observed large-scale distribution 
of the intra-cluster gas (ICG) and its evolution can be directly 
compared to predictions of cosmological structure formation 
models and so constrain them. Also, when using the X-ray/SZ 
method to derive distances to clusters directly, and thus deter- 
mining cosmological parameters, the large scale distribution 
of the ICG has to be known well since incorrect ICG mod- 
els lead to bias in the determinati on of the distance and thus 
in cosmological parameters (e.g., Kawa hara et al.||2008| for 
a summary of sy stematic errors see Molnar, Birkinsh aw &| 
|Mushotzky|2002] ). 

In this paper we focus on what qualitatively new constraints 
on the large scale distribution of the ICG we can expect from 
analyzing data to be taken with the Yuan-Tseh Lee Array 
for Microwave Background Anisotropy (AMiBA; Ho et al. 
2009; Wu et al. 2009) interferometer operating at 94 GHz 
with the planned upgrade to 13 antennas (AMiBA 1 3; Koch et 
al., in preparation). Thus we carry out a feasibility study to es- 
timate how well we should be able to constrain the large-scale 
distribution of the ICG using AMiBA 13. We first summarize 
the presently available observational constraints on the large- 
scale distribution of the ICG (Q. In ^3] we derive a family 
of models for the ICG from our high-resolution cosmologi- 
cal simulations. We then present our methods of generating 
SZ and X-ray images of simulated clusters of galaxies in ^4] 
Our method to simulate visibilities for mock AMiBA13 ob- 
servations is described in ^5]Model fitting and the results are 
presented in ^6] Finally, in^T] we discuss our results for the 
constraints onthe shape parameters of our ICG models from 
mock AMiBA 13 observations. We quote all errors at 68% 
confidence levels (CLs). 

Our companion papers describe the details of the de- 
sign, performance, and the science results from the first ob- 
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servational season of AMiBA with the first configuration 



(AMiBA7). |Ho et al | ( |2QQ9| ) describe the design concepts 
and specifications of the AMiBA telesc ope. Technical as- 
pects of the instrument s are described by |Chen et aL] (|2009 ) 
and |Koch et al.| ([2009 ). Details of t he first SZ obse rv ations 



and data analysis are presented by jWu et al.j ( |2009| ). |Nish 
ioka et al7] ( |2009| ) assess the int egrity of AMiBA 7 data per- 



forming several statistical tests. |Lin et al.| ( |2009j > discuss the 
system performance and verification. Contamination from 
foreground sources and the primary cosmic microwave back- 
ground (CMB) flu ctuations is estimated by Liu et al. (2009). 
[Koch et aLlp010|) present a measurement of the Hubble con- 



stant using AMiBA7 and archival X-ray data. Umetsu et al. 
(2009 ) determine gas mass fractions using gravitational lens 



ing and A MiBA7 observations of galaxy clusters. Hu ang et| 
al. ( [2010 ) discuss cluster scaling relations between AMiBA7 
and X-ray data. 

2. CONSTRAINTS ON THE LARGE-SCALE DISTRIBUTION OF THE 
INTRA-CLUSTER GAS 

Thermal bremsstrahlung is generated by the scattering of 
two particles (an electron and an ion) in the ICG, thus the 
X-ray emission in massive (hot) clusters, where this is the 
dominant emission process, is proportional to the square of 
the electron density, since the ion density is proportional to 
the electron density (e.g., Sarazin 1988). The SZ effect, the 
inverse-Compton scattering of cold photons of the CMB by 
electrons in the hot ICG, is proportional to the electron density 
(Sunyaev & Zel'dovich 1980; for recent reviews see Rephaeli 
1995; Birkinshaw 1999; and Carlstrom, Holder & Reese 
2002). As a consequence, X-ray observations are more sen- 
sitive to the inner parts of clusters, while SZ observations are 
relatively more sensitive to the outer regions. The X-ray flux 
is dominated by signal from regions with 0.1 ^ r/R^ < 0.4, 
where the virial radius, 7? v i r , is defined here according to the 
usage in Bryan & Norman (1998), and the total SZ decrement 
is dominated by signal from regions near 7? v i r (see Figure 10 
ofFang&Haiman2008). 

We measure the SZ signal, which is proportional to the line- 
of-sight (LOS) integral of the electron pressure, and the X-ray 
surface brightness, which is proportional to the LOS integral 
of the X-ray emissivity. The projections cause the observ- 
ables to depend on the LOS size of the cluster. This size can 
be estimated based on the angular size from the X-ray or SZ 
image and the angular diameter distance to the cluster, D A . 
The X-ray and SZ brightnesses of the cluster then provide two 
equations from which we can estimate two unknowns, Da and 
some characteristic electron density in the cluster. Observing 
a sample of clusters, we can derive Da as a function of the red- 
shift, z, and thus constrain cosmological parameters. This is 
usually called the SZ-X-ray (SZX) method (e.g., Birkinshaw 
1999; see Koch et al. 2010 for an application of this method 
using AMiBA7 observations). 

In practical implementations of the SZX method, we deter- 
mine the spatial model for the cluster from the higher signal- 
to-noise (SN) X-ray observations, which typically go out to 
about half of the virial radius. The caveat of this method is 
that there is no guarantee that the ICG distribution at large 
radii follows an extrapolation of the distribution determined 
from X-ray data. Also, due to projection effects, measurement 
errors, etc., the distribution of the ICG determined from X-ray 
measurements might be biased. Models used to describe the 
X-ray observations of clusters are typically (3 models for the 
density distribution, p g ex (1 +r 2 / r core 2 )~ 3/3 / 2 , where the spatial 



extent is determined by the core radius, r core , and the fall off 
by the exponent, (3 (Cavaliere & Fusco-Femiano 1976); with 
either constant temperature (isothermal /3 models), or a grad- 
ually declining temperature as a function of the distance from 
the cluster center. The resulting ft parameters are in the range 
of 0.5 - 0.8 typically. Many relaxed clusters have ft « 2/3, 
which provides a shallow density profile, p ex r~ 2 , at large 
radii (e.g., Sarazin 1988; for recent results see Maughan et 
al. 2008 and references therein). 

However, numerical simulations as well as SZ and X-ray 
observations s uggest a much steeper fa ll off of the density 
at large radii. |Roncarelli et al.| ( |2006| ) used a sample of 9 
clusters of galaxies in the mass range of 1.5 x 10 14 M - 
3.4 x 10 15 M from smoothed particle hydrodynamic simu- 
lations (SPH) to derive gas density and temperature profiles 
in the outskirts of clusters. They used simulations with and 
without cooling, supernova feedback and thermal heat con- 
duction and found that the profiles steepen as a function of 
radius. They also found that cooling and supernova feedback 
do not affect the density and temperature profiles at large radii 
significantly. Their results support the theoretical expectation 
that the distribution of gas at large radii in clusters of gal axies 
is determined mainly by gravity. Hallm an~et al.[ ( [2007] ) fit- 
ted isothermal ft models to mock X-ray and SZ observations 
of simulated clusters drawn from adaptive mesh refinement 
(AMR) cosmological simulations. They found that isother- 
mal ft model fits to X-ray surface brightness distributions of 
simulated clusters are biased to low ft values relative to fits 
to SZ distributions, and that the fitted ft val ues depend on 
the pr ojected outer cut off radii used. When Hallman et al. 
(2007) used a projected radius limit equal to R Y [ Y , the ft pa- 
rameters based on SZ structures scattered around (3=1. Haug- 
bolle, Sommer-Larsen & Pedersen (2007) derived an empiri- 
cal model for the pressure distribution in clusters of galaxies 
based on high-resolution SPH simulations and observations. 
They also found a steeper fall off of the pressure at large radii 
t han that predicted by X -ray observations. 

|Afshordi eta l. (2007) used WMAP 3 year data to stack im- 
ages of 193 massive clusters of galaxies and detected the SZ 
effect s tatistically out to about 2 R Y j Y . Using a larger cluster 
sample, [Atrio- Barandel a et aL] (|2008) determined the average 
electron pressure profile in clusters by stacking 700 known 
clusters extracted from the 3 -year WMAP data. They showed 
that the isothermal (3 model does not p rovi de a good fit on 
large scales . Both |Afshordi et aL ( 2007| ) and Atrio-Barandela 
|et al.| ( [2008] ) concluded that an ICG model with a density pro- 
file with a fall off of p g oc r~ 3 at large radii and a temperature 
profile derived from hydrostatic equilibrium is a good descrip- 
tion of their data. 

The large-scale distribution of the ICG was studied in 
three individual clusters of galaxies (Abell 1835, Abell 1 914 
and CL J1226.9+3332) by Mr oczkowski et al | ([2009) us- 
ing SZA observations at 30 GHz (and at 90 GHFlor CL 
Jl 226.9+3332). The SZA is an i nterferometer consisting 
of eight 3.5 meter diameter dishes (Muchov ej et al.| |2007). 
They used a parameterized pressure profile with five param- 
eters based on cosmological numerical simulations of Nagai 



et al (2007). |Mroczkowski et al.| ( [2009l ) fixed the three slope 
parameters at their values derived from simulations and X- 
ray observations, and fitted only for the amplitude and the 
pressure scale radius using their SZ data. They used a den- 
sity distribution derived from X-ray observations to deter- 
mine the temperature distribution based on the ideal gas law 
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Fig. 1 . — Spherically averaged gas density, p, (in units of the critical density, p c ) and temperature, T (in keV) distributions of the relaxed massive clusters CL1 
(left panels) and CL2 (middle panels) and non-relaxed cluster CL3 (right panels) as a function of radius in units of the virial radius, R ylv (dashed lines). The error 
bars represent rms variations due to spherical averaging. The solid lines show the best-fit density and temperature models to the respective 3D distributions. 



(temperature ex pressure/density). Mrocz kBwski et al.| ([2009) 
found that the SZ profiles drop more steeply than predicted by 
isothermal f3 models, and, similarly to previous studies, that 
the profiles drop more steeply than predicted by a /3 = 2/3 
model even if the change in the temperature is taken into ac 
c ount (Figure 3 of Mroczko wski et al.|2009' 



IVikhlinin et al.| ( 2005| ) measured the temperature profile in 
\i low redshift relaxed clusters using Chandra data. In 
three clusters the temperature profiles were measured out to 
about 0.7# vir . At r > 0.1-0.2 R vir they found that the fall off 
of the temperature with radius is sel f similar in relaxed clus- 
ters when scaled by R Y j Y . Recently Ettori & Balestra (2009) 
and |Bautz et aL] (|20()9 ) studied the out er regions of galaxy 
clusters using X-ray observations. Ettori & Balestra (2009) 
used Chandra observations of 11 clusters with SN ratio 
greater than 2 out to r > 0.77? v i r . The low count rate in the 
outer regions of galaxy clusters did not allow them to deter- 
mine the temperature distribution out to R Y i Y . They derived 
the slope of the gas density and temperature distribution at 
the virial radius assuming hydrostatic equilibrium. Ettori & 
Balestra found that the X-ray surface brightness distribution 
is steepening with larg er radii, implying a n equivalent f3 « 1 
(within errors) at R vir . |Bautz et aLl f2009) used SUZAKU ob- 
servations of relaxed cluster Abell 1795. They mapped the 
X-ray surface brightness and temperature distribution out to 
about 0.9 Rvi Y in two directions and found f3 = 0.64 within 
r < 1 Mpc. At larger radii they found a steeper fall off of the 
density in the South (/? > 0.64), but a rising density profile 
towards the North with a maximum at 1.9 Mpc (1R Y [ Y ). The 
increase of the X-ray surface brightness in the North direction 
might be due to a contribution from a filament in the LOS. 



3. MODELS FOR THE ICG FROM AMR SIMULATIONS 

We derive self- similar spherically symmetric models for the 
distribution of the ICG in relaxed clusters of galaxies using a 
sample of clusters drawn from cosmological AMR simula- 
tions performed with the cosmological code ENZO (O'Shea 
et al. 2004) assuming a spatially flat cold DM model with 



cosmological parameters (Q m , ft\, Qb, h, cr 8 ) = (0.3, 0.7, 
0.047, 0.7, 0.92), where Q m , Q,^ encode the current mat- 
ter and baryon densities and the cosmological constant, ag is 
the power spectrum normalization on 8 fr 1 Mpc scales, and 
the Hubble constant Hq = 100 h km s -1 Mpc -1 . This cosmo- 
logical model is close to the model implied by the WMAP 5 
year resu lts except for cr 8 , which is much larger (Dunkley et 
|al.|2009| ). The AMR simulations were adiabatic (in the sense 
that no heating, cooling, or feedback were included). The box 
size of the original, low resolution, cosmological simulation 
was 300 h _1 Mpc. The clusters of galaxies in our sample were 
re-simulate d with high resolution usin g the same technique as 
described in Younger & Bryan ( 2007]). 

The resolutions (minimum AMR cell size) of the high- 
resolution simulations at R = 0, 1 and 4R Y i Y were about 25 kpc, 
80 kpc and 250 kpc. The total virial masses of the 10 massive 
clusters in our cluster sample fell between 1 and 2 x 10 15 M . 
Relaxed clusters were selected based on their density distribu- 
tion: after the removal of filaments, we chose clusters with a 
smooth spherically averaged density profile with little angular 
variation, no sign of recent major merger events, and a relaxed 
core (for more details about our simulated cluster sample and 
analysis see Molnar et al. 2009). Out of a total of 10 clusters, 
two clusters satisfy our criteria for relaxed clusters (CL1 and 
CL2). We show spherically averaged gas density and temper- 
ature profiles for our massive relaxed clusters and one massive 
cluster with a non-relaxed core, CL3, in Figure [T] We include 
CL3 to check if we could constrain the large scale distribu- 
tion of the ICG in a cluster which has a non-relaxed core, but 
is otherwise relaxed. The error bars represent the rms of the 
density and temperature variations due to spherical averaging. 
While the density distributions are similar in all relaxed clus- 
ters, showing only small deviations from radial averaging, the 
temperature profiles show more variation. This is due to the 
sensitivity of the temperature to shocks from merging and in- 
ternal flows. The solid curves in Figure [T] show the density 
and temperature profiles of the best-fit models for all simu- 
lated clusters (see below). The physical parameters of our 
selected clusters are summarized in Table 1 . 

We use spherically symmetric double fj models truncated at 
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Fig. 2.— Likelihood contours (68%, 95.4% and 99.7% CLs, solid red, 
dashed green, and dash-dotted blue lines) as a function of the outer (3 model 
parameters for fitting double (3 models to the 3 -dimensional density distribu- 
tion of CL1, 2, and CL3 (top to bottom). The square and triangle in the top 
figure and the squares in the middle and bottom panel represent best fit val- 
ues from fitting non-isothermal double (3 models to simulated X-ray images 
in Projection X and Y of CL1, Projection Y of CL2 and Projection X of CL3. 

the virial radius, R ylv , to describe the density distribution of the 
ICG in massive relaxed clusters. We use the same f unctional 
form for the temperature distribution at large radii as Loken et 
aL]([2002 ), and a gaussian to describe the core region. The gas 
aensity and temperature models for massive clusters within 
the virial radius can be summarized as: 



p{r) = p l l {l + r 2 /r 2 l T^ 2 + p 2 2 {l + r 2 /r ( 



.2 r 3/3/2 
core/ 



and 



T(r) = T (a c exp[-r 2 /r c 2 ] + 1) (l + r/r T ) 



(1) 



(2) 



where the large scale distribution is described by r core , /3, rj 
and 5, and a possible extra component at the center is parame- 
terized by r\, fa, a c and r c , p\ and p2 are the central densities, 
and the transition between the two fj models at a \ is facilitated 
by 6\(r,a\) = 0(a\-r) and a{) = Q(r-a\), where is the 
Heaviside step function. We determine the best-fit parameters 
for each cluster by maximizing the likelihood functions 



[(^)-(^m)/}] 2 



where (Fo)i is the median value of the 3-dimensional (3D) 
density or temperature (F = p or T) of a simulated cluster 
in the z-th radial bin, (F M )i is the corresponding value pre- 
dicted by the model considered, and api is the correspond- 
ing standard deviation (here we use an approximation and as- 
sume that the sum of the fluctuations due to substructure, as- 
phericity, etc. are gaussian, which is a reasonable assumption 
since we excluded filaments with large positive density fluc- 
tuations). These functions provide good fits to the density and 
temperature profiles of our selected clusters out to the virial 
radius: the fits are well within the la error bars due to spher- 
ical averaging, except for one point where the deviation is la 
for a temperature profile (see solid lines in Figure [T). The 
best fit parameters are summarized in Table 1. In Figure 
we show the likelihood contours as a function of r core and 
(the shape parameters of the second j3 model, which describe 
the large scale distribution of the ICG) to the 3D distribution 
of the density in the outer parts of our selected clusters (the 
other parameters marginalized). The contour levels were de- 
termined based on the standard \ 2 statistic. The likelihood 
contours show how well constrained the parameters are, sub- 
ject to our assumption of spherical symmetry, and allow us to 
estimate the level of degeneracies between parameters. In car- 
rying out fits to the 3D temperature distribution of simulated 
relaxed clusters we noticed that the exponent of the tempera- 
ture model, 5, does not change much from cluster to cluster (in 
agreeme nt with an analysis using more clusters by |Loken et| 
|al.|2002| ). Also, the best- fit models are not significantly better 
than models with fixed 5= 1.6. Thus we fix 5 = 1.6 in our fit- 
ting and in the rest of our analysis. Since the best fit f3 values 
are close to 1 in all selected clusters, we conclude that ICG 
models with fj = 1 and S = 1.6 provide good fits to these clus- 
ters including CL3, which has a non-relaxed core. This result 
verifies our assumption that the outer region of CL3 is relaxed. 
Therefore, we find that the density distribution at large radii 
can be approximated with a power law, r~ a , with a = 3, which 



■Xf, 



(3) 



is close to a = 3.4 as found by |Roncarelli et al.] (|2006). Our 
results suggest that the pressure (oc ~pT), at large radii can be 
approximated with a power law with a = 3(3 + 5 = 4.6, which 
is close to a = 5, fo und by Na gai et aL] ( |2007| ) and used by 
|Mroczkowski et al.] ( |2009| ). We conclude that the density and 
temperature functions (Equations p] and H\ are adequate for 
relaxed clusters and provide a family of ICG models that can 
be fitted to observational data. 

X-ray observations of clusters of galaxies show that most 
relaxed clusters have cool cores (e.g., Vikhlinin et al. 2006), 
suggesting that cool-core clusters are relaxed. Cosmologi- 
cal numerical simulations suggest that early major mergers 
destroy the developing cool cores, but cool cores are strong 
enough to survive late major mergers, and thus cool cores are 
associated with cluster formation history and not necessarily 
with the dynamical state of clusters (Burns et al. 2008). This 
conclusion seems to be supported by the observational result 
that a substantial number of relaxed clusters do not posses a 
cool core at their center (e.g., HIFLUGS sample^]). In our 
sample of clusters we do not have cool core clusters, but our 
method would work with either cool-core or non-cool-core 
clusters since we model the core of the cluster separately. In 
the case of cool-core clusters the amplitudes of our central 
temperature model, a c , would be negative. Since outside of 
core gravitational physics dominates (Roncarelli et al. 2006), 
and the core region is modeled separately, we conclude that 

8 http://www.astro.virginia.edu/ cls7i/papers/HIFLUGCS_CC.pdf. 
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TABLE 1 

Fitted shape parameters to clusters of galaxies from AMR simulations 



ID a M vir b R yk c ai d n d ft r cor e d P a c r c d n d S e 

CL1 9.1E+14 2.0 0.05 0.013 0.10 0.113 1.033 0.26 0.052 1.07 1.6 

CL2 1.2E+15 2.2 0.08 0.025 0.13 0.139 1.001 0.57 0.011 1.40 1.6 

CL3 f 1.1E+15 2.1 - - - 0.110 1.096 0.46 0.056 1.00 1.6 



a Galaxy Cluster ID 

b Virial mass in Solar Mass 

c Virial radius in Mpc 

d in units of R vlT 

e fixed 

f single p model 

our adiabatic cosmological simulations and cluster models are 
adequate for our purpose. 



4. SZ AND X-RAY IMAGES OF SIMULATED CLUSTERS 

We derive the 2-dimensional SZ and X-ray surface bright- 
ness distributions for simulated clusters in projections along 
the X, Y and Z axes (projections to the YZ, XZ and XY 
planes, Projection X, Y and Z, hereafter). We ignore relativis- 
tic effects, which is a good approximation in our case since the 
intra-cluster gas temperatures in our relaxed AMR clusters is 
less than about 10 keV. We derive the SZ signal in Projection 
Z by integrating along the LOS (I which coincides with z in 
this case) over the extent of the cluster (from £\ to £2) using 



ATc L (x,y) = ATcwNc l LSZ / n e (x,y,£)T e (x,y,£)d£, (4) 

where x and y are spatial coordinates in the plane of the sky, 
i.e. perpendicular to £\ n e = f g p g / n e mp is the electron den- 
sity, where p e is the mean molecular weight per electron 
and m P is the proton mass; p g is the gas density; f g is the 
mass fraction of baryons in the cluster that are contained in 
the ICG (we adopt f g = 0.9), and we use the standard as- 
sumption that the electron temperature equals the gas tem- 
perature, T e = T. The frequency dependence is contained 
in ATclo = p(*v) Tqmb ^bcft / (m e c 2 ), where the dimensionless 
frequency x v = /^/(^b^cmbX where Tqmb is the monopole 
term of the CMB, h P and k B are the constants of Planck and 
Boltzmann, and the function p(x v ) = x v coth(x zy /2) -4 (e.g., 
Birkinshaw 1999). The SZ normalization is 

Nclsz= / n e (0At)T e (0At)d£. (5) 
Jii 

In practice, we pixelize x, y and £, and approximate the inte- 
gral with a sum over the LOS from £\ = -10 Mpc to £2 = 10 
Mpc. Similar expressions were used for Projections X and Y. 

Liu et al. (2010) studied the contamination from CMB, 
galactic diffuse emission and point sources in six clusters of 
galaxies observed in the first year of AMIBA. They found 
that the contamination is dominated by CMB fluctuations. 
The low level of contamination by point sources at around 
90 GHz is due to the low angular resolution of AMIBA and 
the falling spectra of most radio sources. Even though some 
sources have inverted spectra, theoretical predictions based 
on VLA observations at lower frequencies suggest that only 
about 2% of clusters are contaminated at a signif icant level by 
point sources at this frequency (see Figure 13 of Sehgal ^taL] 



|2010| ). Therefore we include CMB contamination in our SZ 
image simulations but not point sources, since we can select 
relaxed clusters with no significant point source contamina- 
tion for structural studies. Thus we have 

AT(x,y) = AT CL {x,y) + ATcMB{x,y), (6) 

where the CMB contamination, ATqmb, is generated as in 
Umetsu et al. (2004). 

We derive the X-ray surface brightness in Projection Z by 
integrating along the LOS {£) over the extent of the cluster 
(from £\ to £2) as 

rh 

Scdx,y) = S C ix>Ncix / n 2 e (x,y,£)T l J 2 (x,y,£)d£, (7) 

where Sew is the central X-ray surface brightness and the nor- 
malization, Nx, is 

Nclx= n](W)T^\W)d£, (8) 

where we use the scaling Sx oc y/T e for thermal 
bremsstrahlung (similar expressions were used for Pro- 
jections X and Y). Similarly to the SZ surface brightness, we 
pixelize x,y and £, and approximate the integral with a sum 
over the LOS from £ x = -10 Mpc to £ 2 = 10 Mpc. 

We simulate X-ray images of our relaxed AMR clusters as- 
suming a field of view (FOV) of 16' x 16' with a pixel size of 
2" x 2" . We sample the Poisson distribution with expectation 
values equal to the integrated flux per pixel as 

Fx(x,y) = F CL (x,y)+B x , (9) 

where we calculate FcL(x,y) using Equation [7] assuming 
Few = 7.5 cnts/pixel (Sew in units of integratea flux/pixel) 
and a uniform background of B x = 0.2 cnts/pixel (typical pa- 
rameters for X-ray observations of clusters of galaxies). 

The Compton-y images, y = j dr (k^T e )/{m e c 2 ), where r 
is the optical depth, the SZ images (cluster plus CMB at the 
AMiBA frequency band, 94 GHz), X-ray surface brightness, 
and simulated X-ray images including background noise, of 
our two relaxed clusters (CL1 and CL2) and one non-relaxed 
cluster (CL3) in Projections X, Y and Z assuming that they 
are located at a redshift of 0.3, are shown in Figures [3] H] and 
[5] The Compton-y, X-ray surface brightness and simulated 
X-ray images are shown in logarithmic scale (first, third and 
fourth rows), and the SZ images with CMB contamination 
(second row) are shown in linear scale. The virial radii of our 
massive relaxed clusters are about 2 Mpc, which span about 8 7 
on the sky at this redshift. The dark blue regions -1 .2 mK) 
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Fig. 3. — Simulated images for Projections X, Y, and Z of relaxed cluster 
CL1 (left, middle and right column). From top to bottom: Compton-_y, AT 
(including the cluster SZ signal and CMB fluctuations at 94 GHz), X-ray 
surface brightness and simulated X-ray image including background noise 
(logarithmic color scale except SZ map which is linear). The image size is 
50 ' x 50 1 for the SZ images (first two rows) and 16' x 16' for the X-ray 
images (last two rows). The physical scale of 1 Mpc, same within each raw, 
is represented by horizontal bars. 



on the SZ images mark the cluster centers, the yellow and red 
regions represent positive and negative CMB fluctuations with 
an amplitude of about ±130/iK and an rms of about 90 /iK. 
In the Compton-y images we can follow the diffuse gas out 
to about 3R Y { Y , w here the external s hocks of massive clus- 
ters are found (see Mol nar et al.|2009| ), but in the SZ images 
the diffuse gas around clusters seems to extend out to about 
R w i Y only, due to contamination from CMB fluctuations, which 
dominate the large-scale structure. However, the SZ images 
of relaxed clusters, in most projections, show similar charac- 
teristics within the virial radius: a circularly- symmetric cen- 
ter and somewhat elongated outer regions. The core regions 
in the X-ray and SZ images (for example regions with yellow 
color in rows 1 and 3), in a few projections (Projection Z of 
CL1 and Projections X and Z of CL2), show asymmetry due 
to asphericity of the cluster and contamination by filaments in 
the LOS. The core region of non-relaxed cluster CL3 seems 
to be disturbed in all projections. 



5. AMiBA13 VISIBILITY SIMULATIONS 

Interferometers measure visibilities, the Fourier transforms 
of the intensity distribution multiplied by the primary beam of 





FIG. 4. — Same as Figure [3] but for simulated images of relaxed cluster 
CL2. 

the telescope. In the small-angle approximation, the visibili- 
ties can be written as 



/CO POO 
/ 
00 J— t 



B v (x,y)I v (x,y)e 



-il7r(ux+vy) 



dxdy, (10) 



where V u (u, v) is the visibility function in the uv plane, which 
is the Fourier conjugate of the positions, x and y on the sky. 
B u (x,y) and I v (x,y) are the primary telescope beam pattern 
and source intensity at x and y at frequency v, and we ignore 
effects due to a finite bandwidth, finite time averaging, and 
other effects associated with practical interferometers. We 
convert temperature differences to intensity units using 



A/ 



^CMB 



X 1/ £ 



AT 



(11) 



CMB 



where Iqmb is the intensity of the monopole term in the CMB. 
The visibilities, Vql and Vcmb at frequency v are calculated 
from Ar C MB and AI C l using Equations [T0| and |TT) 

We simulated visibilities for mock two-patcnAMiBA13 
observations of our relaxed clusters assuming a compact con- 
figuration for the 13 antennas (Figure [6} for a detailed de- 
scription see Koch et al., in preparation). In this configuration 
the many short baselines provide the highest sensitivity to the 
large-scale radio structure. The instrument noise and the con- 
tribution from the CMB were simulated in visibility space and 
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Fig. 5. — Same as Figure[3]but for simulated images of cluster CL3. 

added to the visibilities of the simulated clusters (see Umetsu 
et al. 2004). We assumed an observing time of 60 hours on- 
source. The errors in the azimuthally averaged visibilities for 
the simulated clusters, for the CMB and for the instrument 
noise are all comparable at this exposure time. Using longer 
exposure would reduce the errors for the instrument noise, but 
would not lower the errors either for our ICG models or for 
the CMB. Therefore we would not be able to reduce the er- 
ror bars significantly by using longer observations. Usually 
a subtraction scheme is used in cluster SZ observation with 
half of the time on-source, and half of the time off source to 
reduce systematics, ground pick up and sky background (e.g., 
Wu et al. 2009). Therefore 60 hours of on-source integration 
involve 120 hour of total observing time. This can be conve- 
niently accommodated into a 1 -month observing schedule. 

We use only projections of relaxed clusters which seem to 
show no contamination by filaments or distortion due to pro- 
jection effects (Projections X and Y of CL1 and Projection Y 
of CL2), and a non-relaxed cluster, CL3, in a most relaxed 
projection, projection X. We derive simulated visibilities for 
two-patch observations as 



V(Ui,Vi) = V sour ce(Ui,Vi)-V bkng (Ui,Vi) : 



(12) 



where the on-source visibilities, V sc 



are 



and the off-source (background) visibilities, Vbknd, are 

Vbkndfe, V f ) = Vc MB (u h vO + Vnoisefe, V,-), (14) 

where Vsz is the visibility of the SZ signal, Vcmb, Vnoise* and 
^cmb ' Koise are two sets °f visibilities of the CMB fluctuations 
and noise (different for on- and off-source observations). We 
assumed that the CMB fluctuation fields for the two patches 
(about one degree apart) are uncorrelated, which is a conser- 
vative assumption: this way we somewhat overestimate the 
noise due to the CMB. Since our models are spherically sym- 
metric, we have no constraints on them from the imaginary 
part of the model visibilities (the imaginary part is identically 
zero), therefore we work only with the real part of the visi- 
bilities. Note that in real applications the imaginary parts can 
be used to check the amplitude of the CMB fluctuations in 
the field, assess non- sphericity of the cluster and the pointing 
accuracy. 

As an illustration, in Figure [7] we show the radial profile of 
the real part, TZe(V), of the azimuthally averaged simulated 
visibilities of one realization of a two-patch observation for 
CL1 in Projection Y for the 90 and 98 GHz AMiBA chan- 
nels (Channels A and B, Ho et al. 2009) as a function of 
R uv . Since the length scale for the visibility data is in units of 
the observing wavelength, the visibilities are represented by 
two sets of curves. The red diamonds and green squares with 
error bars represent visibilities for Channel A and B. The er- 
ror bars represent the rms of the azimuthally averaged real 
part of the visibilities at the AMiBA13 baselines. Visibilities 
of other projections of CL1 and projections of CL2 and CL3 
with no contamination from filaments are very similar due to 
the structural similarity among the cluster SZ images in these 
projections (see Figures [3j [4] and [5]). 

6. MODEL FITTING 

We use non-isothermal double f3 models truncated at the 
virial radius to describe the ICG. We determine the best- 
fit parameters using likelihood functions. Our model for 
the SZ surface brightness distributions (Equations [T] and 
[2]) has eleven free parameters: ten shape parameters, p = 
(fli,ri,^i,r core ,^,fl c ,r c ,rT,^,^ V ir) 5 and one amplitude, ATq. 
Unfortunately, due to limited spatial resolution, FOV and 
receiver noise, we do not expect to be able to determine 
all ten parameters using AMiBA 13. Thus, similarly to 
Mroczko wski et aL] ( |2009| ), we reduce the number of free pa- 
rameters m our models. We proceed the following way: first 
we determine the shape parameters for the central part of the 
cluster, a\,r\, f3\,a c ,r c , and the core radius for the large scale 
distribution, r core from simulated X-ray observations (fixing 
P = 1, r T = 1.0 Ry{ Y and 6 = 1.6 based on our results from nu- 
merical simulations). We assume that the X-ray surface den- 
sity and temperature are determined with acceptable accuracy 
only in the central part of the cluster, to 0.57? v i r and 0.27? v i r . 
Therefore r T cannot be determined from fits to the X-ray tem- 
perature profile, and the X-ray emissivity within 0.57? v i r gives 
no useful constraints on the large scale distribution of the tem- 
perature. Having determined the shape parameters for the 
central part, we derive constraints on the temperature scale 
parameter, r T , from simulated AMiBA 13 visibilities. 

We calculate the predicted integrated X-ray flux/pixel for 
our X-ray model as 



Vsource(Wi,Vf) = VciX^, Vf) + VcMBte, Vf) + V n oise(^, Vf), (13) 



FxM - FdM + BxM , 



(15) 
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AMiBA dish configuration: N diah = 13 
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x/cm 

Fig. 6. — Antenna layout for AMiBA13 in the compact configuration (13 
dishes with 1.2 m diameter) used for simulations in this paper. Other possible 
antenna positions are marked with x-s. 



where F C lm is our cluster model and B XM is the background. 
We calculate Fqlm as 



F C lm(x, y, P) = FiN X ilxi(x, y, P) + F 2 N X2 Zxi (x,y,p), (16) 

where F\ and F 2 are the central integrated flux/pixel for 
the two model components, the normalizations are = 
IxiQO ,0, p) and A^ 1 =X Z2 (0,0,p), and 



X xl (x,y,p) = 2 / F 2 {l + r/r T y 5 l 2 dl 



Pic 

X X2 (x,y,p) = 2 / T 2 (l + r/r T y^ 2 d£, (17) 
Jo 

where r 2 = x 2 +y 2 +f, and the cut off in the LOS is £ 2 = 7? 2 ir - 
x 2 -y 2 . We ignore the error in the X-ray background, B X m 
since it is negligible relative to other sources of error, such as, 
for example, errors due to azimuthal averaging. 

We experimented with the likelihood ratio for Poisson dis- 
tribution and Mighell's statistic (see Mighell 1999 for a 
detailed analysis) in fitting this structure. We found that the 
fitted parameters obtained by these two methods are virtually 
identical. Therefore we decided to use the Poisson likelihood 
ratio test: 

-lnC x = ^Mi-Ni+Ni HNi/Mi), (18) 



where Ni and M t are the observed and expected numbers of 
photons. 

We derive spectroscopic-like X-ray temperature profiles, 
T sp , for our simulated clusters using the weighting scheme of 
Mazzotta et al (2004), which has been shown to provide tem- 
perature profiles similar to those observed with CHANDRA 
and XMM (Nagai et al. 2007), 



/ w sp TdV 



Tsp ~ Jw sp dV 



(19) 
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Fig. 7. — Real part, He {V} (in Jy), of the azimuthally averaged simulated 
visibilities as a function of uv radius (R uv ) in Projection Y of the CL1 (plus 
CMB and noise, one realization) for AMiBA13 Channel A and B (red di- 
amonds and green squares) in a compact configuration (see Figure|6j. The 
error bars represent instrumental errors of 60 hour observations. Tne best- 
fit p model is also shown for the two channels (red solid and green dashed 
lines). 

where the weight is w sp = n 2 /T 3 ^ 4 . We determine the best fit 
parameters for the temperature model maximizing the likeli- 
hood function 



-21n£ r = 



l(T sp )i-(T M )i}? 



(20) 



where (T sp )i and (T M )i are the median values of the observed 
and model spectroscopic temperatures, and a t is the standard 
deviation in the ith radial bin. 

Since the shape parameters of the inner part of clusters, r\, 
Pi, r core , a c , and r c , are not sensitive to the large-scale distri- 
bution of the temperature (heavily weighted towards the cen- 
ter of the cluster), we determine them fixing fj = 1, r T = 1 
and S = 1.6 (based on our results from fitting to the density 
and temperature distributions of our simulated clusters). In 
our case, since (3=1, the X-ray emissivity, which is propor- 
tional to (1 + r 2 /r core 2 ) -3 in the outer parts of the cluster, drops 
about six orders of magnitude from the cluster center to the 
virial radius. Therefore a moderate change in R Y { Y (say 20%) 
causes only an insignificant change in the X-ray signal (Equa- 
tion [17]) due to a change in the upper limit, t c (except around 
R Y { r where the X-ray signal is negligible; see Figures[3]|4]and 
[5]). For our purposes, therefore, we fix the value ofR y { r as- 
suming that an estimate for its value with a 20% accuracy is 
available from other measurements. We determine the shape 
parameters, r\, f3\, r core , a c , and r c by maximizing the like- 



lihood functions, Equations 18 and 20 In practice, since the 



spectroscopic-like temperature is not sensitive to a few tens of 
percent change in the /3 model parameters, we can determine 
the shape parameters using iteration. We use the likelihood 



function for the X-ray emission (Equation 18) to determine 
r\, /3\ and r core , and then the likelihood function for the X-ray 
temperature (Equation [20]) to determine a c , and r C9 then we 
iterate over these steps (usually only two steps are needed). 
This method proved to be faster than a search for the maxi- 
mum likelihood in the five dimensional parameter space. 

In Figure [8] we show the radial profiles of the X-ray surface 
brightness distribution, the spectroscopic temperature and the 
SZ signal for Projections X and Y of CL1, Projection Y of 
CL2, and Projection X of CL3 (blue plus signs with error bars 
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Fig. 8. — Radial profiles of X-ray surface brightness, gas temperature, normalized SZ surface brightness and gas density (top to bottom) in projections X and Y 
of CL1, projection Y of CL2 and projection X of CL3 (left to right columns). Blue plus signs and error bars: AMR cluster in projection (no noise, no CMB; first, 
second and 3rd rows); red squares: simulated data points for one realization including contamination; black plus signs and error bars: spherically averaged gas 
temperature and density distributions in simulated clusters based on their 3D distributions (shown out to the virial radius; second and 4th row); blue dashed lines: 
best fitted models to simulated X-ray and AMiBA13 observations (first, second and 3rd rows); blue solid lines: 3D distribution of the temperature and density in 
simulated clusters based on fits to the projected (2D) X-ray surface brightness, gas temperature, and SZ surface brightness profiles (second and 4th rows). 



representing dispersion in azimuthal averaging). For each 
projection, we also show one realization of our Monte Carlo 
simulations (red squares), and, for comparison, the tempera- 
ture and density profiles derived from their respective 3D dis- 
tributions (black plus signs with error bars representing dis- 
persion in spherical averaging, solid black lines represent best 
fit models; same as in Figure [TJ. We also show the best fit X- 
ray surface brightness models (blue dashed lines) derived as 
described in the previous paragraph. The other dashed and 
solid blue lines, explained later in this section, are included 
here for later convenience. Based on the first row in this fig- 
ure, we conclude that our non-isothermal double f3 models 
provide good fits to the X-ray surface brightness profiles in 
our relaxed clusters (CL1 and CL2), and even in our cluster 
with a non-relaxed core (CL3). 

The best fit core radii of the outer f3 model from fitting our 
non-isothermal double f3 models to the X-ray surface bright- 
ness distributions and temperature profiles (with /3 = 1 fixed) 
in different projections are shown in Figure [2] (squares and 
triangles). These values are within the 68% CL of the best fit 
values obtained from fitting double (j models directly to the 
3D distribution of the density (solid red lines). We conclude 
that the core radius for the large scale distribution can be de- 
termined accurately, with no bias, from X-ray data using our 



ICG models. 

As a second and final step, we calculate the visibilities for 
our models, V M , at frequency v, using Equations [T0| and [TT 

as 



Ar CLM (x,y) = AT N sz lsz(x,y,p), 



(21) 



where ATq is the central SZ amplitude, the normalization is 

N^=lsz(0,0,p) 9 md 



P1+P2 



(l + r/r T y 5 dt (22) 



Note, that since r/r core and r/r T are both dimensionless, 
Ar CLM (x,)0 depends on R W { Y only through i c . We use Equa- 
tion [12] to derive simulated observed visibilities, Vo, for each 
realization as described in Section ^5] 

We determine the best fit parameters for our models by 
maximizing the SZ likelihood function defined as 
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Fig. 9.— Likelihood contours (68%, 95.4% and 99.7% CLs, solid red, 
dashed green, and dash-dotted blue lines) for fitting for the normalized SZ 
amplitude (ATq/ATq^) and the temperature scale radius rj in projections X 
and Y of CL1, projection Y of CL2 and projection X of CL3 (from top to 
bottom). 



-2lnCsz = J2 



[Ke{V (R l uv ; vj)} - 1Ze{V M (R l uv ; ^)}] 2 



4 



(23) 



where 1Ze{Vo} and lZe{Vu} are the azimuthally averaged real 
part of the observed and model visibilities at the ith uv ra- 
dius, R l uv , at frequency Vj, where j = 1,2 for the two AMiBA 
frequency channels, and cr^- is the Gaussian noise, which is 
assumed to be the same for all antenna pairs and frequencies. 

Similarly to the X-ray signal, a moderate change in R Y [ r 
(all the other parameters fixed) causes only an insignificant 
change in the SZ signal (Equation 22) because the pressure 
is about four orders of magnitude smaller at the virial ra- 
dius than at the center of the cluster. In this final step, we 
assumed that we determined the model parameters for the 
central region of clusters (n, /3\, a c , r c ) and r core from fit- 
ting our models to the X-ray surface brightness and tem- 
perature distribution (as above), therefore, among the shape 
parameters of our non-isothermal double /? models, p = 
(n , P\ , ^core , 1 , «c , r c , rx , 1 .6, R Y i Y ), there is only one unknown 
parameter, rj. We determine the best fit for r T and the normal- 
ized SZ amplitude (A7^ z /A7q) by maximizing the SZ likeli- 
hood function, Equation [23 



In Figure M we show the 2D projected temperature and SZ 
profiles for tne best fits (dashed blue lines, 2nd and 3rd row). 
Solid blue lines represent the 3D best fit radial profiles for the 
temperature and density (2nd and 4th row). In some cases the 
deprojected temperature profiles are underestimated by 10- 
20% close to the cluster center because cold substructures in 
the LOS reduce the spectroscopic temperature and the weight- 
ing suggested by Mazzotta et al (2004) does not always cor- 
rect for this effect accurately. However, these deviations in 
the central temperature do not affect the fits to the large scale 
SZ profiles as observed by AMiBA 1 3 because of the low res- 
olution of the instrument. 

In Figure [9] we show the CLs for the temperature scale ra- 
dius r T and the normalized SZ amplitude in different projec- 
tions of X and Y of CL1, Projection Y of CL2 and Projec- 
tion X of CL3. The normalization, ATq/ A7o£, is the ratio of 
the best-fit Ar for each Monte Carlo simulation to the best- 
fit value determined from all Monte Carlo simulations. The 
confidence levels have been determined using Monte Carlo 
simulations. The contours for the 68% CL, for example, are 
the smoothed version of the contours containing 68% of the 
simulated best fit points based on simulations. 



7. DISCUSSION 

We have simulated AMiBA13 observations of massive re- 
laxed clusters of galaxies including CMB contamination and 
receiver noise using clusters drawn from cosmological numer- 
ical simulations to assess how well we should be able to con- 
strain the large-scale distribution of the ICG. Our simulated 
SZ images (rows 2 in Figures [3] [4] and [5]) show that at 94 GHz, 
on the scale of the cluster cores (few arc minutes), the clus- 
ter SZ signal dominates the CMB fluctuations, so that CMB 
contamination is not important in surveys searching for rich 
clusters. On a scale of ten arc minutes, corresponding to the 
extent of the ICG in massive clusters at z ~ 0.3, the CMB 
contamination is at a comparable level to the SZ signal. Con- 
tamination from CMB fluctuations is also important in the re- 
gions in galaxy clusters close to their virial radius. Spectral 
separation of the cluster SZ and the CMB signals based on 
multi-frequency observations seems to be essential for study- 
ing the outskirts of galaxy clusters and the SZ signature of the 
large-scale structure. 

Using our AMR simulations, we showed that a spherical 
non-isothermal double f3 model with a temperature distribu- 
tion described by Equation [2] provides good fits to the radial 
distributions of the ICG in our selected massive clusters. We 
generated X-ray and SZ images of our clusters drawn from 
numerical simulations assuming that the clusters are at a red- 
shift of 0.3. 

We used the simulated X-ray data to determine the shape 
parameters, r\, f3\, r core , of the double fj model for the den- 
sity distribution and the central model parameters, a c and r c of 
the temperature model by minimizing the likelihood functions 
Equations [18] and [201 We assumed that R yiv is constrained 



to ±20% by other measurements, for example, from gravi- 
tational lensing, as done by Broadhurst and Barkana (2008) 
and Umetsu et al. (2009). This is a conservative estimate 
since combining weak and strong l ensing data, one can deter- 
mine R yiv with about 5% accuracy ( Umetsu et al.|2010| ). The 
values of the likelihood functions corresponding to the 68%, 
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95.4% and 99.7% CLs for fitting for r core and /3 to the 3D den- 
sity distribution of simulated clusters are shown in Figure [2] 
The elongated shape of the confidence levels in this figure 
shows a degeneracy between the two shape parameters of the 
f3 model: r core and f3. Combinations of small (r core ,/3) or large 
( r core,P) both give good fits to the density distributions. A 
similar degeneracy has been reported for fitting /? models to 
X-ray observations (e.g.,|Grego et al.|2001[ |Reese et aL) 2000 ; 
Birkinshaw & Hughes|1994| ). 

As a final step, we used simulated AMiBA13 visibilities to 
determine the temperature scale radius, r T , and the SZ nor- 
malization (fixing all other parameters) by maximizing the 
SZ likelihood function, Csz (Equation [23]). Our results are 
shown in Figure [9] For model fits to Projections X and Y 
of CL1, we obtain r T = 1.023^4 #vir, and r T = 1.60^;§ # vir . 
Fits to Projection Y of CL2 yield r T = 0.927^ #vir. We ob- 
tain r T = 0.631^o 2 ^vir for our cluster with a non-relaxed core, 
CL3, in Projection X. 

All results for the temperature scale radius, rj, are within 
68% of the best-fit values based on fitting the 3D distributions 
(Table [T}. The 68% CLs for all clusters with circular pro- 
jected X-ray distribution (Projection X of CL1 and Projection 
Y of CL2) and the cluster with a non-relaxed core (Projec- 
tion X of CL3) are within 50% of the best fit values. There- 
fore our results suggest that, using relaxed clusters with cir- 
cular morphology, we should be able to use AMiBAl 3 along 
with X-ray observations to obtain unbiased parameters for our 
non-isothermal double (3 models (even for clusters with a non- 
relaxed core, such as CL3), and constrain r T within 50%. Our 
model fits to Projection Y of CL1, a relaxed cluster with ellip- 
tical morphology, although returning r T within 68% of the 3D 
fitted value, shows a large error due to asphericity. We expect 
that using a more accurate elliptical model would result in bet- 
ter constraints on rj. The SZ amplitude, which is important 
for the determination of the Hubble constant for example, can 
be determined with 3-4%, which is better than the expected 
accuracy of the absolute calibration of AMiBAl 3. 

We carried out simulations at different redshifts between 
0.1 and 0.4. For a distant cluster the beam dilution reduces 
the signal; if the cluster is too close, the outer parts of the 
cluster fall outside of the FOV. We have found that the optimal 
redshift for determining the large-scale distribution of the ICG 
with AMiBAl 3 is z ~ 0.3. 

We conclude that we should be able to use AMiBAl 3 to 
determine the large-scale distribution of the ICG in massive 
relaxed clusters of galaxies located at a redshift of 0.3 by de- 
termining the temperature scale radius with an about 50% sta- 
tistical accuracy. AMiBA, as upgraded to 13 dishes with 1.2- 
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m in diameter, will be a powerful tool for constraining the 
large-scale distribution of the ICG. 

The degeneracy between density and X-ray temperature can 
be broken and determined out to the virial radius using cur- 
rent X-ray telescopes. However, an accurate determination of 
the temperature profile out to requires a long exposure 
time (due to the low photon count rates at the outer parts of 
clusters) and mosaicing (due to the limited FOV). This is the 
reason why only a few attempts have been carried out to map 
the outer regions of clusters (see Constraints on the large 
scale distribution of the ICG similar to AMIBA can be derived 
using bolometer cameras. In principle, bolometer cameras 
with arc minute resolution can cover clusters with the nec- 
essary sensitivity out to R yir (ACT: Hinks et al. 2009; SPT: 
Plagge et al. 2009), although this has not yet been done. 

In this paper we considered spherically symmetric galaxy 
cluster models, and focused on statistical errors due to the 
AMiBAl 3 telescope and receiver system. We should be able 
to reduce the observation time needed to measure rj by using 
a more sophisticated observational strategy. A single point- 
ing to the center of the cluster, although it simplifies the data 
analysis, has the disadvantage that it has a reduced sensitivity 
at the outer parts of the cluster, where the signal is weaker. 
We expect that mosaic observations including pointings to- 
wards the outer regions in clusters, although more difficult to 
analyze, would enable us to reach our goal with an exposure 
time shorter than 60 hrs. We leave simulations to quantify the 
effects of mosaic observations and more sophisticated ICG 
models on parameter determination for future work, as well 
as a detailed study of parameter estimation from a large set 
of simulated relaxed clusters with a wider range of mass and 
redshift. 
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